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Abstract 

Recent experiments on short MoGe nanowires show a sharp superconductor-insulator transition tuned by the normal 
state resistance of the wire, with a critical resistance of R c m Rq = /i/(4e 2 ). These results are at odds with a broad 
range of theoretical work on Josephson-like systems that predicts a smooth transition, tuned by the value of the 
resistance that shunts the junction. We develop a self-consistent renormalization group treatment of interacting phase- 
slips and their dual counterparts, correlated cooper pair tunneling, beyond the dilute approximation. This analysis 
leads to a very sharp transition with a critical resistance of Rq. The addition of the quasi-particles' resistance at finite 
temperature leads to a quantitative agreement with the experimental results. This self-consistent renormalization 
group method should also be applicable to other physical systems that can be mapped onto similar sine-Gordon 
models, in the previously inaccessible intermediate-coupling regime. 



Key words: 

PACS: 74.78.Na, 74.20, 



z, 74.40.+k, 73.21. Hb 



1. Introduction 

One of the most intriguing problems in low- 
dimensional superconductivity is the understanding 
of the mechanism that drives the superconductor- 
insulator transition (SIT). Experiments conducted 
on quasi-one-dimensional (ID) systems have shown 
that varying the resistivity and dimensions of thin 
metallic wires can suppress superconductivity [112] , 
and in certain cases lead to an insulating-like be- 
havior [3141516171 . 

Particularly interesting are recent experiments 
conducted on short MoGe nanowires [6] that explore 
the SIT tuned by the wire's normal state resistance 
with a critical resistance R c « Rq . Resistance mea- 
surements of the quasi-lD MoGe nanowires reveal 
a strong temperature dependence that can be fitted 
with a modified LAMH theory [819] of thermally ac- 
tivated phase-slips down to very low temperatures. 



However, for these narrow wires it appears that the 
LAMH theory is valid only in a narrow temperature 
window (see discussion in Sec. [5]). Moreover, the 
LAMH analysis does not explain the appearance of 
a critical value of R c f=a Rq . 

The universal critical resistance may suggest 
that, at a temperature much lower than the mean- 
field transition temperature, T <C T c , the wire acts 
as a superconducting (SC) weak link resembling a 
Josephson junction (JJ) connecting two SC leads. 
Schmid [TU] and Chakravarty [TT] showed that quan- 
tum phase-slip fluctuations in such a JJ lead to a 
SIT as a function of the junction's shunt resistance, 
i? s , with a critical resistance of Rq, and that the 
resistance across the junction obeys the power law 



R(T) cx T V " a J . The theory was later extended 
to J J arrays and SC wires [12113114] , Within these 
theories, a similar power-law prevails. However, 
contrary to this general prediction, Bollinger et 



Preprint submitted to Elsevier 



2 February 2008 



al. [B] observe that the resistance of the MoGe wires 
exhibits a much stronger temperature dependence, 
even close to the SIT. 

In a previous work |15j , we have presented an ap- 
proach that captures both the critical resistance of 
R c « Rq at the SIT and the sharp decay of the 
resistance as a function of temperature. We treat 
the SIT in nanowires as a transition governed by 
quantum phase-slip (QPS) proliferation. This pic- 
ture alone, however, cannot account for the observed 
strong temperature dependence of the resistance. 
We argued that the key ingredient left out in pre- 
vious works is the inclusion of interactions between 
QPSs in such a finite-size wire, especially when the 
phase-slip population is dense. 

We treat these interactions in a mean-field type 
approximation: when analyzing the behavior of a 
small segment of the wire, we include in its effec- 
tive shunt-resistance the resistance due to phase- 
slips elsewhere in the wire. This scheme is motivated 
by numerical analysis of a related problem, an in- 
teracting pair of resistively-shunted JJs [16]. This 
self-consistent treatment primarily produces a sharp 
temperature dependence of the resistance. In addi- 
tion, we include the effects of the Bogoliubov quasi- 
particles, which couple to the potential gradient cre- 
ated by each phase-slip [17] ■ Consequently, the re- 
sistance obtained in the experiment can be fitted 
without resorting to the LAMH theory beyond its 
limit of validity. 

In this manuscript we generalize our previous re- 
sults to the weak Josephson coupling limit, where 
conductance through the wire proceeds by Cooper 
pair tunneling. The use of a self-consistent treatment 
of correlated Cooper pair tunneling events leads to 
a similarly sharp SIT in the limit of highly resistive 
wires. This self-consistent approximation of phase- 
slip interactions should be applicable to similar mul- 
tiple sine-Gordon models in the theoretically chal- 
lenging intermediate-coupling regime. 

The remainder of the manuscript is organized as 
follows. In Sec.[2]we discuss the validity of the LAMH 
theory for the wires in Ref. [6j. In Sec. [3] we sum- 
marize our previous results on the role of quantum 
phase-slip interactions in short wires, with details of 
derivation of the microscopic model in Appendix |A"1 
We generalize these results to highly resistive wires 
in Sec. [4] and present our conclusion in Sec. El 



2. On the validity of the theory of thermally 
activated phase-slips in short wires 

In an attempt to explain the strong temperature 
dependence of the resistance of quasi-lD MoGe 
nanowires, Bollinger et al. [BJ have shown that the 
experimental curves can be fitted with a modified 
LAMH theory [8|9] of thermally activated phase- 
slips, down to very low temperatures. Nevertheless, 
estimations based on the parameters of the wires 
in Ref. [6] (Table QJ suggest that most of the tem- 
perature range in the experiment lies outside of the 
region in which the theory is applicable. The the- 
ory of thermally activated phase-slips is based on 
the time-dependent Ginzburg-Landau description 
of a superconducting (SC) wire. This description is 
valid at temperatures higher than the gap, and far 
enough from T c , such that fluctuation corrections 
are small, T* < T < T G . Here T* is defined by 
A(T*) = T*, with A(T) the temperature depen- 
dent order parameter, and Tq = T c (l — Gi), with 
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the Ginzburg-Levanyuk num- 
ber for the quasi-lD wires, where = Rw£,/L 
is the normal resistance of a section of the wire of 
length £. For the wires in Ref. [6J, the LAMH the- 
ory is valid only in a narrow temperature window 
as T* w 0.9T C , and estimates for Tq range between 
0.82T C and 0.97T C (see Fig. [Hand Tabled]) [15115] . 
Moreover, the LAMH analysis does not explain the 
appearance of a critical resistance R c « Rq . 

This estimate suggests that most of the relevant 
temperature range in the experiment conducted in 
Ref. [6] is in the regime T <C A(T), where quantum 
fluctuations become increasingly important. While 
the universal critical resistance supports the idea 
that the transition occurs due to quantum phase 
slip proliferation in a narrow constriction, the naive 
extension of the theory to SC wires [12113] cannot 
account for the observed sharp drop of the resis- 
tance. In the following section we will present a self- 
consistent treatment of quantum phase-slip prolif- 
eration, that captures both the critical resistance of 
Rc ~ Rq at the SIT and the sharp decay of the re- 
sistance as a function of temperature. 



3. Low-resistance wires 



The microscopic action for a SC wire can be ob- 
tained from the BCS Hamiltonian by a Hubbard- 
Stratonovich transformation followed by an expan- 
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Fig. 1. Validity of the LAMH theory for the data of wire 
No. 6 of Ref. [6] (see Table [T] for detailed parameters of the 
wire). The fitting parameters used to fit the data are the 
coherence length, § = 6.7 nm, and the critical temperature, 
T c = 4.7 K, marked by an arrow. These fitting parameters 
yield a temperature range T* = 4.23K < T < T G = 4.55 K, 
for which the LAMH theory is expected to be valid, shown 
in gray (see the text for a definition of T* and Tq). 

sion around the saddle point (see Appendix [20 . 
In the limit of low energy scales, u>, Dq 2 <C Ao, this 
yields 

L l/T 

S = N Q AA 2 Jdx J drj^- [ln(p 2 ) - l] 
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where L and ^4 are the wire's length and cross sec- 
tion, respectively, = ttD/8A ,v p — \J (37r/2)£>A 
the amplitude velocity, — ^/ ttDAq(2 AV c Na + 1) 
the phase velocity, V c the Fourier transform of the 
short-range Coulomb interaction, Nq the density 
of states, D the electronic diffusive constant in the 
normal state, and the SC order parameter is pa- 
rameterized as A = Ao/oe 1 ^, with Ao the mean- field 
solution. For the wires in Ref. [6], 2AV c Nq oc N± ~ 
1000 ^> 1, leading to v p -C v$ oc v p \/Nj_. Here 
■Wj- = PfA/tt 2 is the number of ID channels in the 
wire. 

This action supports QPS excitations, which are 
characterized by two distinct length scales: v p /Aq oc 
£ <C £,y/N± oc v<f,/Ao. For very long wires, £ <C 
£-\/W± C i, in the dilute phase-slip approxima- 
tion, this problem can be mapped onto the pertur- 
bative limit of the 1 + 1-dimensional sine-Gordon 
model 12 J. In the opposite limit of very short wires, 



L < £ <C £y/N±, the system resembles a J J and can 
be mapped onto the + 1-dimensional sine-Gordon 
model. However the wires in Ref. [6] appear to be in 
the intermediate regime, £ -C £ -C £\/N_\_. Hence, 
while phase-slips occur in different sections of the 
wire, they are indistinguishable, as each creates a 
phase fluctuation that spreads over distances larger 
than the wire itself. 

Moreover, the wires in Ref. [6] have a sizable bare 
fugacity. Using Eq. |T]), one can estimate the core 
action of a phase-slip of duration tq = 1/ Ao- Choos- 
ing the following trial function for a phase slip: 



<t> = arctan (v^t/x) 

P = min [\J (x/x ) 2 + (A t) 2 , 1] , 



(2) 



and identifying the part of the action [Eq. (JTJ], 
that corresponds to p ^ 1 as the core action, we 
minimize this expression with respect to the phase- 
slip diameter, xq. Keeping only leading terms in 
N± = (pj^A/n 2 ), this yields a core action of S c — 
71-/8^/1/3(1 - e)R Q /R i 22J. Here e = 1 - T/T c is 
the reduced temperature and = R-w^/L is the 
normal resistance of a section of length £. Using 
this expression, the measured values of Ryj and L 
(Table [l}, and the assumption that £ ss 20 nm, the 
phase-slip fugacity is estimated to be ( — e~ 5c s» 
0.05 — 0.45 for the different wires. Consequently, 
in the critical region, i?vv Rq, there is a dense 
population of phase-slips that interact with one an- 
other; thus, the dilute phase-slip approximation is 
no longer a proper description. 

The flow equation for the fugacity of a phase-slip 
anywhere in the wire is given by: 



dl 



Rq 
R* 



C, 



(3) 



where dl = —din A, and A is the running RG scale. 
Eq. Q treats the phase-slip as occurring on an ef- 
fective J J, with R s being the effective shunting re- 
sistance of the entire wire. When £ is small, R s will 
include only the effective impedance of the leads. If 
£ is not very small, we will need to include in Eq. ([3]) 
additional terms of higher powers of £, which de- 
scribe interactions between QPSs. To deal with a fi- 
nite £, we include the resistance due to other phase- 
slips in the wire in the effective shunt resistance of 
the junction, i? s [T5]. This is akin to guessing the 
form of a complete resummation of higher-order £ 
terms in Eq. © [23|24j . 
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This treatment was successfully tested by numer- 
ical analysis of a simpler analog of the system, an 
interacting pair of resistively shunted JJs [16) . In 
this work, Werner et. al studied the phase diagram 
and critical properties of the pair of JJs using Monte 
Carlo simulations and renormalization group calcu- 
lations. The authors found that, in the region of the 
intermediate coupling fixed point, there is a remark- 
able resemblance in the critical behavior between 
the two-junction system and a single junction. In or- 
der to explain this resemblance, Werner et. al sug- 
gest that the two-junction system can be described 
by an approximate mean-field theory. In the mean- 
field approximation, each junction at criticality be- 
haves as an independent junction and sees the other 
junction as an effective resistor whose resistance is 
determined by phase-slip events. This picture can 
account for the observed (and calculated) properties 
of the two junction-system at the superconducting 
to normal transition point. 

The main physical intricacy of the self-consistent 
approach is the determination of the effective shunt 
resistance, R S (Q), that governs the renormalization 
of the phase-slip fugacity [Eq. (J3j)] . A phase-slip pro- 
duces time- varying phase gradients, and hence elec- 
trical fields. These dissipate through two channels in 
parallel: the SC channel - which has an effective re- 
sistance due to other phase slips, R ps - and the quasi- 
particles conduction channel [17j . which has resis- 
tance i? qp . Once the disturbance reaches the leads, it 
also dissipates through the electro-dynamical modes 
of the large electrodes, whose real impedance is pa- 
rameterized by i? e lec- 

For T <C T c , the resistance of the quasi-particles, 
i? qp , can be approximated by 



m L n 
qp = ~ T = n 

e T n Tiqp /i ^qp 



T 



2ttA 



(4) 



Unfortunately, we lack a microscopic model for 
the impedance of the electrodes, as this depends 
on the details of the system such as the junction's 
shape and material. However, we expect that at 
large scales, T < A < A, the electrodes will act as a 
transmission line to the electromagnetic waves gen- 
erated by the phase-slip. This transmission line is 
characterized by a real impedance which we denote 
as i? c iccj and use as a fitting parameter. Hence, the 
effective shunt resistance that affects the renormal- 
ization of the phase-slip fugacity at T < A < A [Eq. 
©] is 



R S [C(A)]=R ( 



lec 



1 



i 



i? ps [C(A)] i? qp (T) / 



(5) 



Fig. [5] shows the circuit we suggest describes the 
system. 
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Fig. 2. The effective electrical circuit for the nanowire. A 
phase-slip produces time- varying phase gradients, and hence 
electrical fields, that dissipate through the quasi-particle con- 
duction channel, which has resistance -R qp , and through the 
SC channel, that has an effective resistance due to other 
phase-slips, -R ps . Once the disturbance reaches the edges, 
it also dissipates through the electro-dynamical modes of 
the large electrodes, whose real impedance is represented by 

^clcc ■ 

The resistance is measured in response to an ap- 
plied DC current. In this zero frequency limit, the 
electrodes act as a capacitor connected in parallel 
to the wire. Therefore, the measured resistance is 
the total wire resistance, unaffected by the environ- 
ment, which is cut off from the wire: 



Rtot(T) = (l/iwcco] + i/flqpCnr 



(6) 



The occurrence of a phase-slip causes a resistance 
in the otherwise SC wire through the relation R ps tx 
(L /£)C 2 - Using this relation and Eq. ©, we can 
write an RG equation for the dimensionless resis- 
tance 



dl * dl 



2 1 



R 



Q 



d(R ps /R Q ) 
dl 



2 1 



Rs(0 

R s (R ps )J \R Q 



Rq 



(7) 



Integration of Eq. ((7]) , with the effective resistance 
given in Eq. from the ultraviolet (UV) cutoff 
A(T*) = T* to the infrared cutoff T, yields R ps {T). 

The wire's DC resistance Eq. ©, calculated using 
Eqs. ([5]) and ([7]), is plotted in Fig[3](a) as a function 
of temperature for different R* . Here R* = R to t {T*) 
is the normal state resistance of the wire, at the UV 
cutoff A(T*) = T* . We assume that the wires are 
thin enough such that the mean-field transition from 
normal to SC is wide, and R* w R to t(T+) = R u . For 
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simplicity, we have assumed throughout our calcula- 
tion that A(T) w Ao. This assumption holds in the 
low-temperature regime, T <C A, where the theory 
of QPSs, based on the effective action Eq. ([J), is ex- 
pected to be valid. The resistance of the environment 
is taken to be i? e iec = 0.1i?Q. In practice, A, R* and 
i?ei ec can be used as fitting parameters. Moreover, 
R* can be determined independently as the resis- 
tance measured below the drop that indicates pass- 
ing through T c of the SC films (see Ref. ©). 




Fig. 3. The total wire resistance, Rtot, in units of Rq 
as a function of renormalized temperature T/A for differ- 
ent R* /Rq, ranging from 0.8 (lower plot) to 1.01 (upper 
plot) in increasing steps of 0.03. The environment resis- 
tance was taken to be R e \ ec /RQ = 0.1. Gray traces compare 
the total wire resistance -Rtot(T) = ^l^?^S P ^J2,s to a sin- 

gle Josephson junction with R(T) = R* (T/A) 2(H « / ' Rs " 1) , 
and constant R s . (a) The total wire resistance, 

ps-Rqp(^) 



Rtot(T) 



r , in the limit of weak phase-slip fugac- 



-Rps+-Rqp(T) ' 

ity calculated using Eqs. J5} and JTJ. The insulating plots are 
cutoff at temperature, To, for which f (To) = 1 [see the discus- 
sion prior and subsequent to Eq. J5)]. (b) The total wire re- 
sistance, Rtot(T) = l/G tot (T) = (G cp (T) + l/R qp (T))- 1 , 
in the weak Josephson coupling limit calculated using Eqs. 
(TTTJ and llOt . Here we have assumed the wire length is 
L = 5£, where § is the coherence length. The initial drop 
in the resistance of the high-resistance wires is a manifes- 
tation of the difference between the total wire resistance 
and the effective shunt resistance, in the presence of local 
quasi-particle relaxation, in the weak coupling limit. As the 
ratio between the wire length and the coherence length in- 
creases, the initial drop in the resistance becomes more pro- 
nounced. The low-resistance plots are cutoff at To, for which 
J (To) = 1. Inset: The resistance of the wire calculated in 
the self consistent approximation, Eq. lilt , shows an initial 
sharp increase as the temperature is lowered, followed by a 
moderate increase. This moderation of the diverging resis- 
tance is due to the finite density of quasi-particles, present 
at finite temperature. 

Eq. |(7J) is also applicable to a wire with R* > Rq. 
In this limit the fugacity increases in the renormal- 
ization process, and Eq. ([3]) is no longer valid for 



£ > 1. For wires with R* > Rq, we overestimate 
the phase-slip fugacity as C(T*) — 0.5 and integrate 
Eq.©, with 



R B [t(A)] = R< 



lec 



+ 



J/ C (fl*)C(A) 2 Rqp(T) 



(8) 



from the UV cutoff, A(T*) = T* , down to the tem- 
perature To for which C(^o) = 1- Here the propor- 
tionality constant y^(R*) is set by the initial condi- 
tion R* = (l/[j, c (i?*)C(T*) 2 ] + l/£ qp (T))-\ with 
C(T*) = 0.5. The results are shown in Fig. [3J (a). 
Overestimating C(T*) gives an upper bound on To, 
where Eq. © is no longer applicable. Fig. [3J (a) 
shows that the transition between SC and insulating 
wires occurs for a critical resistance R c w Rq. How- 
ever, in contrast to the standard Josephson junction 
theory (gray curves in Fig. [3J (a)), the transition is 
much sharper (notice the logarithmic scale). 





Measured Values 


Fitting Values 


Curve 


L(nm 


) Rw(kn) 


R*(k{l) A(K) R 


dec (kfl) 


1 


177 


5.46 


4.2 


2.5 


1.2 


2 


43 


3.62 


2.62 


2.35 


1.25 


3 


63 


2.78 


2.13 


3.07 


0.66 


4 


93 


3.59 


2.89 


3.85 


0.55 


5 


187 


4.29 


4.5 


6.55 


0.31 


6 


99 


2.39 


2.09 


4.84 


0.4 



Table 1 

Summary of nanowire parameters, and the parameters used 
to fit the experimental data; L is the length of the wire 
determined from SEM images, iiw is the wire's normal state 
resistance, taken as the resistance measured below the film 
transition. Fitting parameters: R* = Rtot(T = A) is the 
wire's resistance at the UV cutoff, A the SC order parameter, 
and -R c i ec the impedance of the electrodes at T < A < A. 

A comparison between the theoretical curves and 
the experimental data taken from Ref. [5] is shown 
in Fig. |4j The curves were calculated by fitting R* , 
A, and i? c i oc - Since the theory of QPSs is expected 
to be valid at T -C A, deviations from the theo- 
retical curves at high temperature are reasonable. 
In general, as A is proportional to T c , increasing A 
shifts the sharp decay of the resistance to high tem- 
peratures. Both R* and i? ioc affect the high tem- 
perature resistance, whereas R* and A control the 
width of the transition. 

We have made an attempt to fit the data corre- 
sponding to the insulating wires of Ref. [5] . As the 
insulating wires are thinner, we expect a strong sup- 
pression of T c [2] , which sets the scale for the high 
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temperature cutoff T <C A(T). Moreover, in these 
wires, the bare fugacity is estimated by Q — e~ Sa w 
0.32 — 0.84, which results in a relatively high To, 
for which C(To) = 1. Consequently, the theory of 
QPSs with £ < 1 is valid in a narrow range of tem- 
peratures. While we manage to fit the experimen- 
tal curves in this regime of parameters, we do not 
present the results as our fits cover a range of ~ f 
data points, with three fitting parameters. 
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To use this relation and Eq. ©, we need to con- 
sider what the effective shunting resistance of an 
individual Josephson junction (a segment of length 
£) is. When a Cooper pair tunnels a distance £, a 
dipole of strength 2e£ forms. It can relax either by 
locally fusing back through the quasi-particle chan- 
nel, i.e, through a resistance R qp (T)£/L, or by the 
tunneling Cooper pair, and the hole it leaves be- 
hind increasing their separation, until they leave 
the system through the electrodes. The latter pro- 
cess implies the 2e charge goes through a resistance 

-Rcioc + (g cp [J(A)] + R . J (T) ) . These two chan- 
nels appear in parallel, see Fig. [5] Therefore the full 
shunting resistance is: 



i? s [J(A)]- 1 = 



' n ' eleC + V«qp(' r )Gcp[J(A)] + l J 



^qp(r)f 



(10) 



T(K) 



Fig. 4. Comparison between the theory (solid line) and the 
experimental data [6]. Details of the wires, and fitting pa- 
rameters are summarized in Table[T] Note that the estimated 
values for R* are in good agreement with the measured val- 
ues of the resistance of the wire at T = A. The discrepancies 
between theory and experiment occur for the longer wires 
(curves 1 and 5). In these wires the condition v^/Aq 3> L 
might not hold, and one should consider additional renor- 
malization at higher scales. Inset: The effect of an adjust- 
ment of the fit parameters, in roughly 10% (from dark gray 
to light gray): A = 4.356K , R clcc = 440fi and R* = 188M. 



4. High-resistance wires 

To better describe the insulating wires, we turn 
to the weak Josephson coupling limit where phase 
coherence across the wire is lost, and conductance 
proceeds by means of Cooper pair tunneling events 
across regions of fluctuating order parameter ampli- 
tude. The flow equation for the Josephson coupling 
anywhere in the wire is given by: 



dJ 



=^ = 1 



3l 

Rq 



J. 



(9) 



Once more we expect that the Josephson tunnel- 
ing between two such sections will be effected by 
all other tunneling events in the wire, via capacitive 
coupling. 

The pair tunneling event leads to a finite con- 
ductance through the insulating wire, Gcp J - 
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Fig. 5. The effective electrical circuit for the nanowire in 
the weak coupling limit. When a Cooper pair tunnels a dis- 
tance £, a dipole of strength 2e£ forms. It can relax either by 
locally fusing back through the quasi-particle channel, i.e, 
through a resistance R qp (T)£/ L, or by increasing the sep- 
aration between the tunneling Cooper pair and the hole it 
leaves behind, until they leave the system through the elec- 
trodes. The advancing Cooper pair can relax by means of 
successive tunneling events with resistance 1/G cp , or by ex- 
citing quasi- particles along the wire, with resistance R qp (T). 
Once more, as the excitation reaches the edges, it dissipates 
through the electro-dynamical modes of the large electrodes, 
whose real impedance is represented by -R e lec- 

Thus the RG equation for the conductance, in 
units of Gq = 4e 2 //i, is 





= 2 















(11) 



where R s is given by Eq. (|I0[) . The total re- 
sistance of the wire, R to t(T) = I/G?tot(T) = 
(G cp [J(T)] + l/R qp (T)y\ calculated using Eqs. 
(fTU|) and (fTP]) , is plotted in Fig[3](b) as a function of 
temperature for different R*. As mentioned above, 
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R* = R to t(T*) w R n is the normal state resistance 
of the wire, at the UV cutoff A(T*) = T*, which 
can be determined from resistance measurements 
at a temperature below the SC transition of the 2D 
films [6J . The order parameter is assumed to be con- 
stant throughout the calculation and the resistance 
of the environment is taken to be -R i cc = O.IRq. In 
Fig. [3] (b) we have assumed L = 5£. In Ref. [6] the 
high-resistance wires are typically thinner, and are 
therefore expected to have a strong suppression of 
T c [2] • This will lead to a relatively large coherence 
length, and a smaller ratio L/£. 

The initial drop in the resistance of the high- 
resistance wires is a manifestation of the difference 
between the total wire resistance and the effective 
shunt resistance, in the weak coupling limit. The lo- 
cal quasi-particle relaxation reduces the shunt resis- 
tance below Rq at T ~ A(T) , causing an initial drop 
in the resistance. At lower temperature, the den- 
sity of quasi-particles decreases exponentially, and 
R qp £,/L ceases to be smaller than Rq. This leads to 
a decrease in the number of tunneling events causing 
the total resistance to increase. As the ratio between 
the length of the wire and the coherence length in- 
creases, the initial shunt resistance is smaller and 
the initial drop in the resistance becomes more pro- 
nounced. 

The resistance of the wire, calculated in the self- 
consistent approximation, Eq. shows an initial 
sharp increase as the temperature is lowered, fol- 
lowed by a moderate increase. As the pair tunneling 
drops rapidly to zero at a finite temperature, con- 
ductance persists due to the finite density of quasi 
particles. This residual conductance appears as a 
moderation of the diverging resistance, as shown in 
Fig [3] (b) . Conversely, upon neglecting the interac- 
tions between different sections of the wire, the con- 
ductance due to independent Josephson tunneling 
shows a power law decay as a function of decreasing 
temperature, G cp (T) = 1/R* (T/A)^ RQ ^ Ra ~ 1 \ 
Such a power law behavior leads to a finite conduc- 
tance at finite temperature, which shunts the highly 
resistive contribution of the quasi-particles. 

Eq. (TlTj) is also applicable to a wire with R* < 
Rq. In this limit the pair tunneling increases in the 
renormalization process, and Eq. (fTTjl is no longer 
valid for J > 1. We estimate bare Josephson cou- 
pling as J(T*) = 0.5 and integrate Eq. ((9]) with 



R s (A)=R clcc +lyj(R*)J(A) 2 + ——) ,(12) 



from the UV cutoff, A(T*) = T*, down to the tem- 
perature To for which J(Tq) = 1. Once more, the 
proportionality constant is set by the initial condi- 
tion R* = [yj(R*)J(T*f + l/i? qp (T)] _1 . The re- 
sults are shown in Fig. [3] (b). 

When we use the resistive wire theory, we must 
consider the following caveat: If the Josephson cou- 
pling between segments of the wire is only pertur- 
bative, then the phase of the order parameter fluc- 
tuates locally as a function of time. If the frequency 
scale of this fluctuation is comparable to the order- 
parameter magnitude, then there should also be low 
lying sub-gap density of states in the wire, that we 
don't take into account. The theory in this section 
assumes that the phase fluctuations are sufficiently 
slow such that the sub-gap states can be ignored. 

5. Conclusion 

In conclusion, we studied the effect of interactions 
between QPSs in short SC wires, beyond the di- 
lute phase-slip approximation. Our analysis shows 
that treating these interactions in a self-consistent 
manner produces a sharp superconductor-insulator 
transition with a critical resistance R c ~ Rq, in 
agreement with recent experiments [6j. Moreover, 
we have shown that adding the resistance of the 
BdG quasi-particles leads to a quantitative agree- 
ment with the experimental curves. In the dual 
weak Josephson coupling limit, this self-consistent 
RG treatment produces a similar sharp insulator- 
superconductor transition. The sharp drop in the 
conductance of the insulating wires in this limit is 
shown to be accompanied by a residual conductance 
due to the density of quasi-particles at this finite 
temperature. Our method should be applicable to a 
wider range of physical problems which involve the 
proliferation of topological defects with a sizable 
bare fugacity. In particular, it could be applied to 
the study of a Luttinger liquid with an extended 
impurity [25] . 

We thank E. Demler, and P. Werner. Special 
thanks to A. Bezryadin for making his data avail- 
able to us. This study was supported by a DIP 
grant and by an ISF grant. 

Appendix A. Derivation of microscopic 
action 

Consider a system of electrons in the diffusive 
limit that interact via Coulomb repulsion and the 
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phonon mediated BCS interaction. This is described 
by the Hamiltonian H — Hq + Hint + Hbcs ■ 



a.cr' 

Hbcs = 



|r — r' 



(r) 



where an d are electron creation and 

annihilation operators, the Fourier transform of the 
Coulomb interaction in ID V c (q) ~ logg is taken to 
be constant, and u(r — r,-) is the impurity potential 
at point r due to an impurity at point r^. We assume 
that the impurity potential is S correlated, (u(r)) = 
0, and (u(r)u(r 1 )) = 7v lf S(r — r'), where r is the 
impurity scattering time, and Nq is the 3D density 
of states. Following Ref. [26], we apply a Hubbard- 
Stratonovich transformation in order to rewrite the 
BCS interaction and the Coulomb interaction. The 
partition function becomes 

Z = J DADA*DpD4'D^e~ s , 
where 



S = S 



dx 



|AQr) 
A 



+ - I dxdx' p(x)V c -\x - x')p(x') 



+ i J dxp(x) {yj\(x)^(x) + ^[(x^^x)} , (A.2) 

and x = (r, r). Integrating over the fermionic fields, 
the effective action in Nambu-Gorkov spinor nota- 
tion reads 

S = J dx\- l \A{x)\ 2 -TrlnG- 1 



+- / dxdx' p{x)V c T 1 {x - x')p{x') 



(A.3) 



with 
G~ 1 -- 




and £ = — V 2 /2m — p. One can treat the pres- 
ence of nonmagnetic impurities by including a self- 
energy diagram that describes the dressing of the 



electron line by impurities. As a result, the frequen- 
cies and the order parameter are replaced by u = 
t]lu and A = r)A, with r/ = (l + 1/(2t^lu 2 + Af)J 

|27j . Moreover, when calculating polarization bub- 
bles (see section rA.2[) . one must sum over the im- 
purity ladder (namely impurity lines connecting the 
two Green's functions in the polarization bubble) . 

A.l. Uniform order parameter 

In the case of a uniform order parameter, the ac- 
tion in Eq. (|A.3|) can be greatly simplified. In this 
limit, the order parameter A can be chosen to be 
real and the action reduces to [28] 



S[A] = J dx\- x \A(x 

- T J2' ^ d3p 



(2ny 



In f <~ 2 



£ 2 + A 2 . (A.5) 



A variation of the action in Eq. (|A.3[) with respect 
to A* yields the BCS gap-equation 



A ^-y^TA 2 ' 



(A.6) 



where the sum indicates that the frequencies 
are cut off at the Debye frequency, u>d- In order to 
evaluate the sum in Eq. (|A.5[) . we replace S[A] by 



1A 

S[A] - S[0] = J 



dS[A] 
dA 



(A.7) 



where 

as* [A] 



2A 



, 2ATV' / d 3 p-. 

dA A F (^,2 + f2 



£ 2 + A 2 



We have substracted S[0] to avoid divergences. As 
S[0] is independent of A, this choice does not affect 
our final results. Integrating with respect to A re- 
sults in the following action 



(A.4) S[A] - S[0] : 



A 2 
A 

A 2 
A 



47rA T^ (J 



N A 2 



w=0 
1 

2 



A 2 - 



ln 



2uj d 



(A.9) 



which can be written as 
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S[A] - S[0] 



In 



A 2 

A 2 



- 1 



(A.10) 



A. 2. Fluctuations around the mean-field solution 

The action describing phase and amplitude fluc- 
tuations is obtained by expanding Eq. (|A.3[) around 
the mean-field saddle point solution A — Ao . Divid- 
ing the fluctuations into real and imaginary parts, 
5 A = Al + i At, which are connected to amplitude 
and phase variations, the fluctuations around the 
mean field are given by 



S = I dx 



|A | 



Af + A T 



+- J dxdx' p{x)V~ l {x - x')p(x') - Tr In G^ 1 
-TrJ2 — {GoSG- 1 )" . (A.ll) 



Here 

Go : 

5G- 1 




—iuj + £ 



(A.12) 



Keeping only leading terms in p, A^ and Ay, we 
find that the partition function can be written as 

Z = J DA L DA T D P e- s *t^ AL > AT > p \ (A.13) 



with the effective action 



S eS = ~ J <*z(a l A t p)v~ X {<l) 



(a\ 



An 



, (A.14) 



where we have introduced a shorthand notation 
f dq = ATJ2 U J dq/2ir and A is the wire's cross 
section. The screened potentials are given by 



V-\q) 

4 + 









-- + Ua t a 2 



-n A 



Tp 





IIa-tp 

n 



i 



2AV r 



.(A.15) 



/</>/ 



In the above, Tl a f3 is the polarization bubble ob- 
tained by integrating out the electronic degrees of 
freedom, with the vertices a and /3 corresponding to 
the incoming and outgoing bosonic fields. 

In the dirty limit ql, fir <C 1, where q and fl are 
the transferred momentum and frequency, respec- 
tively, and I and r are the mean free path and im- 
purity scattering time, we find 



II A£At (g,Q) = 7riV T^ 

1 



1 



A 2 "I 



WW 



W + W + Dq' 



n Ar A T (g,^)=7r^ T^ 
1 



1 



A 21 
^0 



WW 



W + W + Dq 2 

n pp ( q ,n) = No - ttNqTJ2 

a? 

i 



WW 



W + W + 
n ATP (q,n) = -nN Tj2 



Dq 2 

A n 



i 



WW W + W 



= -n pAT (g,fi). 



-Dq 2 
(A.16) 



A 2 , 



Here J = u+Q, W = y/u 2 + A 2 , W = y/i 
and D — tv 2 f /3 is the diffusion constant. In the low 
temperature limit, for q <C £ ~ £o an d f2 <C Ao, 
where £o is the zero temperature coherence length, 
these may be approximated by 



/ fl 2 it Da 2 
H AiAi (^)-A-^-^(l + T ^ ! + -^ 

U ArrAT (q,n) - X' 1 ^ -No 



n' 2 

4A 2 



vr Dq 2 
4^o" 



n ATP (g, ft) 



in 2 it Dq 2 . 



2A7 l ±_ 6 A 2 8 A 



Note that close to T c , namely in the limit A <C 
ft, Dq 2 <C T, the polarization bubble Ha l a l = 

II AtAt « -JV (f ~ ln (^f^)) . which 

reproduces the time-dependent Ginzburg-Landau 
(TDGL). The term ~ \ft\ describes dissipation. Its 
emergence is a result of the fact that the polariza- 
tion bubbles describe the response of the electronic 
system, in equilibrium. Hence, we have assumed the 



existence of a relaxation mechanism that allows the 
electrons to return to equilibrium. This assumption 
should be taken with caution when studying trans- 
port close to the phase transition, as typical time 
scales close to the phase transition diverge. 

We perform the Gaussian integration over the 
field p to obtain an effective action for At- This 
leads to an extra term in the At propagator, so that 



1 

A 



(HatpY 



-No - 



TrDq 2 



(2k; 

fl 2 



n 



4 A AA 2 2AV C N Q + 1 



(A.18) 



Finally, the effective action describing real and imag- 
inary parts of the fluctuation in the order parameter 
is 



N / dq 1 + 



<r 2 



irDq 2 



12A 2 



|AL(g)| 5 



<? 2 



A 

\A T (q)\ 2 [. (A.l!)l 



/nDq 2 
V 4A ' 4A 2 2AV C N + 1 



We note that the general expression for the uniform 
fluctuations of the magnitude of the order parame- 
ter, Eq. (|A.10[) . reduces to the form given by Ha l a l 
in Eq. (|A.19|) , under the substitution A — A + A L . 
One can show |30j that the real and imaginary parts 
of the order parameter, A^/Ao and Ap/Ao, are re- 
lated under gauge transformation to the amplitude, 
p and phase, <p of the order parameter, respectively. 
Collecting all terms and writing the action in terms 
of the amplitude and phase, we have 

L l/T 

S = N Q AA 2 [ dx I drl^ [In (p 2 ) - l] 









+e 


'p' 2 + 4" 









(A.20) 



with Vcf, and v p as given in Sec. [3] 
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